Coherent quantum evolution via reservoir driven holonomy 
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We show that in the hmit of strongly interacting environment a system initially prepared in a 
decoherence-free subspace (DFS) coherently evolves in time, adiabatically following the changes of 
the DFS. If the reservoir cyclicly evolves in time, the DFS states acquire an holonomy. 
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Geometric phases and Holonomies are generated by 
the presence of a curvature (the Berry curvature) in the 
Hilbert space. The natural way to gain information on 
this curvature is to perform interferometry on a state 
vector which has been dragged along closed loops in the 
Hilbert space. To achieve this task the methods ex- 
ploited are usually divided into two main categories: in 
the first one, a state is driven around the space by means 
of a coherent time evolution 1]. For example, this can 
be obtained by an adiabatic and cyclic motion, which 
forces the state to explore eigenspaces of a parameterized 
Hamiltonian. The second category consists of all those 
motions that result from non-trace preserving quantum 
operations. For example, a series of von Neumann mea- 
surements can be used to evolve wave functions by pro- 
jecting the state vector onto a sequence of overlapping 
subspaces. As an effect of this measurement-induced mo- 
tion, a quantum system may indeed acquire a geometric 
(Abelian or not) phase The Pancharatnam geometric 
phase is a typical example of this category 0, Q . 

Even irreversible quantum processes such as systems 
interacting with Markovian reservoirs can be used to gen- 
erate geometric evolutions. Many environmental models 
have the characteristic of not affecting some particular 
quantum states when these are lying in suitably "pro- 
tected" subspaces, so-called decoherence free subspaces 
(DFS) 4]. States lying in these subspaces are stationar?/, 
i.e. they do not evolve in time. Having control over the 
reservoir (e.g. through engineered reservoirs) may imply 
an indirect control on the protected states of the sys- 
tem 0. In particular, modifications to the parameters 
of the reservoir may result in a controlled time evolu- 
tion of the protected subspace as a whole. If this motion 
is accomplished in a sufficiently smooth way, it can be 
shown that states lying in this subspace evolve coher- 
ently, thereby acquiring information about the geometry 
of the space explored. 

In this letter we will explore the possibility of gen- 
erating (Abelian and non- Abelian) holonomies due to a 
smooth motion of a DFS in the Hilbert space. It is worth 
stressing that the existence of a time-dependent DFS is 
by no means trivial, and its presence often reflects a sym- 
metry preserving evolution. Under a suitable "adiabatic 
condition" it can be shown that a state lying in a DFS 



remains inside the subspace and, hence, is rigidly trans- 
ported around the Hilbert space together with the DFS. 
The evolution experienced is, in fact, coherent, although 
entirely controlled via an incoherent phenomenon. When 
the DFS is eventually brought back to its initial configu- 
ration, the net effect is an holonomic transformation on 
the states lying in this subspace. Notice that the evolu- 
tion obtained is purely geometrical: no dynamics affects 
the system inside the DFS, hence, no dynamical phase is 
accumulated during the evolution. 

Moreover, we will show that in the first order correc- 
tion to the "adiabatic" evolution the DFS is affected by a 
decoherence process. This effect is unavoidable, at least 
conceptually, as far as a non-trivial holonomy is to be 
obtained. On the other hand, the decoherence time in- 
side the DFS can be made arbitrarily large, in the limit 
of strong decoherence and/or slow evolution of the sub- 
space. In fact, the decoherence affecting the DFS be- 
comes weaker as the decoherence acting on the outside 
world gets stronger. This counter- intuitive effect, is the 
essence of the adiabatic approximation. It is, in fact, the 
fast evolving decoherence process affecting the outside 
world which decouples the latter from the DFS. Alter- 
natively, this evolution can be understood in terms of a 
"Zeno effect," where the action of strong environment 
can be regarded as a measuring apparatus continuously 
monitoring the DFS of the system. By using this effect 
in degenerate subspaces it is possible to realize robust 
holonomic gates for quantum computation, analogous to 
those realized by means of adiabatic evolutions ^] . 

Let us first consider a system described by the density 
operator p evolving under the effect of a Markovian envi- 
ronment. The decoherence process due to the interaction 
with such an environment is described by the following 
master equation (h = 1): 



dp 
'dt 



-i[H,p]-J2{^l^kP + prlTk -2Tkprl}, (1) 



fe=l 



where the commutator generates the coherent part of the 
evolution and the remaining part represents the effect of 
the reservoir on the dynamics of the system. The ac- 
tion of each F^ (the Lindblad operators) amounts for the 
different decohering processes that can affect the system. 
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In the Markovian formalism, a decoherence free subspace 
(DFS) is defined as the common eigenspace of all the 
Lindblad operators: Span{|^) | Tk\i{}) = Ck\ip)yk}. 

Suppose that the evolution is solely due to the action 
of the environment, i.e. the presence of a Hamiltonian 
evolution can be neglected compared to the effect of the 
environment {H = 0). Consider, now the situation in 
which the environment depends on some time-dependent 
external parameter. This time dependence might be re- 
garded as either the action of an experimenter who can 
control some degree of freedom of the reservoir, or "drift" 
of some parameter intrinsically occurring in the environ- 
ment. In general, we will assume the rate of this evolution 
very small compared to characteristic time scales of the 
system. Under this assumption, the evolution is governed 
by a time-dependent master equation of the form: 

^ = -^{Ti{t)T,{t)p + prUt)T,{t)-2T,{t)prl{t)}, 

k 

(2) 

where the operator Tk{t) are time dependent. For each 
instant of time we can consider the eigespaces IC{t) — 
SpanjlV') I rj.(t)|'0) = Ck\ip),yk} of the instantaneous 
operators Tkit). Let Il{t) be the projector onto IC{t), 
and nj_ = 1 — n(t) the projector onto the orthogonal 
complement ICj_(t). Consider the operator 

Dit) = ^ {T,it)^T,it) - 2clit)T,it) + \c,it)\H) , (3) 
fc 

whose Hermitian part P{t) = {D{t) + D''{t))/2 is positive 
semidefinite. The space IC{t) can be defined as the kernel 
of D{t), i.e. D{t)Il{t) = 0. The assumptions that we will 
consider are the analogous to the hypotesis of the usual 
adiabatic theorem 0,0|: (i) Il{t) is a chain of projection 
operators, smoothly depending on time, (more precisely 
n(t) norm- twice differenciable 10|) ; (ii) the non-zero 
eigenvalues of P{t) are bounded from below by a time 
independent value 7 > 0. 

Under these assumptions, we will show that: (a) it is 
possible to formulate an adiabatic limit, for which JC is 
decoupled from the environmental action, and therefore 
it is decoherence free; (b) in the adiabatic limit, if the 
subspace K,{ti) at some time ti ^ to + T coincides with 
/C(to) at the initial time to, the neat evolution experi- 
enced by a state belonging to JC is purely geometrical, i.e. 
the evolution between to and ti is described by the holon- 
omy (Abelian or not) associated with the path traversed 
by 1C{t) in the Hilbert space; (c) the first order correction 
to the adiabatic limit is a source of decoherence, whose 
time scale is inversely proportional to the smallest time 
scale of the decoherence affecting the JC±^ , the latter being 
the orthogonal complement of /C. 

We would like to describe the evolution at some time 
ti = to+T of states initially prepared inside the subspace 
)C{to), in the limit in which the time rate of motion of 
)C{t) is much smaller than the smallest characteristic time 



scale of the system, which is 7" . It is then convenient to 
introduce a time independent parameter s = {t — to)/T, 
with s G [0, 1] and solve the master equation in the limit 
of T — > 00. By following a similar approach of the stan- 
dard adiabatic theorem (9|, let's consider a unitary oper- 
ator 0{s) for which: ^ (O^'IIO) = with 0(0) = 1. The 
operator fulfilling this condition is defined by the equa- 
tion z40(s) = G{s)0{s) (with initial condition O^0) = 
1) dll^l where the generator is G = + Q{s), 

and Q{s) is an arbitrary block diagonal Hermitian oper- 
ator, i.e. n(3(s)n_L = 0. By definition, describes the 
change of picture to the " rotating frame" , which rigidly 
follows the subspace /C(s). In fact, in this picture, the 
projector n(s) = 0'< {s)uls)0^sy = n(0), and the cor- 
responding subspace K is time independent. The free- 
dom in the choice of the generator G(s) corresponds to 
the arbitrary choice of a smooth chain of basis within /C 
and IC±, and it is often referred in this context as gauge 
freedom. Under re-parametrization t s — {t ~ to)/T, 
and change of frame p ^ p = pO, equation ^ takes 
the form: 

^ = *[G, p] - ^T{D{s)p + pD^s) - 2 ^ f fcpf t } (4) 

fc 

with G ^ O'^GO, f fe = {O'^VkO - cfcl) /V7, and D = 
O^DO/^. Notice that Ffcll = and that the operator D 
has been re-normalized, so that its Hermitian part P — 
{D + D^)/2 — f ^Ffc has minimum no- null eigenvalue 
not smaller than 1. 

From equation Q it is possible to see that the incoher- 
ent part of the evolution afi^ects only system states lying 
outside the subspace IC. However, this does not prevent 
a state initially prepared in /C to be exposed to deco- 
herence. Indeed, the off diagonal terms Gog = Hj^GII 
and G^jj of G can couple the subspace K. to its comple- 
mentary, and may eventually spoil the coherence of the 
former. In order to unwind the effect of the perturbation 
introduced by G from the main decoherence process, it 
is convenient to introduce a transformation which allows 
the evolution to be described in terms of two effectively 
decoupled manifolds. To this end, consider the following 
effective non-Hermitian Hamiltonian: 

H{s) = -e*^(") {r]G{s) + iD{s)) e"'^'-'^ (5) 

where 77 = l/"fT is the adiabatic parameter and S{s) is a 
(non-Hermitian) operator defined by the condition that: 

(a) pHn_i =_nj_i?n = o (with n_L = i - n) and (/?) 

HS'n = H^S'n^ = [13. Let's now analyze the time 
evolution of the density matrix under the transformation 
e^^ . To this end consider the operator p — e^^ pe~''^\ 
Notice that although this transformation is not unitary, 
p is still a positive semi-definite operator, and hence it 
is a valid density matrix, up to a normalization factor 
N = Tr{pe^^^~^ ^) > 0. By assuming for simplicity that 
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the time derivative of S is neghgible, the evolution of p 
is given by: 

with Tk = e^'^FfcC"*'^. The advantage of this expression 
is that H is now block diagonal and the main coupling 
effects are only due to the second term. We are interested 
in the evolution of the system in the limit of 77 <gC 1. It is, 
then, convenient to consider the expansion of S in series 
of 77: S = T]Si + rfS-z + 0(77'^), and, correspondingly, the 
expansion of the effective Hamiltonian H = Hq + rjHi + 
rf'H2 + 0{rf ). By using the above conditions (a) and 
(/3) it is possible to show that: 



Hi 



-iD{s), 

-nGn-n_LGni, 



H2 = -\n[iSi,G]ii-\ni_[iSi,G]ii^, 

where Si = D-'^GoS - H.c. Notice that 



(7) 
(8) 
(9) 

the 



last expression is well defined, as it is restricted to the 
subspace /C^ where D is invertible. In the limit of ry <C 1 
the evolution of p can be expressed as follows: 



dp 
ds 



1 



C-i[p\ + Co[p\ + vCi[p\ + o{-n' 



(10) 



The largest contribution in H10() is given by 



k 



{vlXuP + pVlXk 



2ffcpfi}, (11) 



where we also expressed f ^ = f ^ +ir][Si,fk\ + 0{rf) and 
retained only the zero-th contributions in rj. By definition 
of f eq. Hll|l acts trivially on the partial density operator 
PDF = npll, i.e. C-i[pdf] = 0. Therefore, the most 
relevant non-trivial evolution affecting the state poF is 
given by the term £0 in the adiabatic expansion: 

Cq[pdf] = i[Gdf,Pdf] with Gdf = HGn, (12) 

where the fact that TkpDF — has been used. This 
expression explicitly demonstrates that retaining only 
the terms up to the zeroth order in the adiabatic pa- 
rameter rj yields to an evolution for poF which is uni- 
tary^ and, therefore, coherent and trace preserving. This 
shows the existence of an adiabatic limit for 77 — > 0, 
in which the evolution is confined - by the decoherence 
process itself - into a suitable subspace, in which the 
evolution maintains its coherence. It is worth notic- 
ing that, due to the non-unitary transformation e*"^, 
the new decoherence-free subspace K, is slightly modi- 
fied from the original K,. Indeed, the presence of the 
generator of the adiabatic motion, G, deformes the DPS, 
whose (generalized) projectors are to be identified with 



n ^ n 



JS 



He 



n -|- i77{5'i, n}. As expected, in 



the limit r/ — > 0, /C converges to K.. 



The evolution of the partial density matrix pjj p can be 
easily solved in the adiabatic limit and formally expressed 
as: 



p^U{s)p{Q)U\s) with U{s) =:Pexp 



« / GDF{T)dT 



where 7^ is a path ordering operator. Let's now consider a 
closed evolution of /C, i.e. an evolution for which n(ti) = 
n(to)- Then it is possible to express the total evolution 
as: 



U{1) ^Vexp 



I I Gdf 



•{T)dl 



^Vcxp (j> AdX 



where, the right hand side expresses the time evolution 
in terms of a path integral. The operator A, is a (non- 
Abelian) holonomic connection, defined as a vector of 
components 



Ax^ ^ -no 



fdO{Xi,\2 



-n, 



(13) 



where A's are a set of variables parameterizing the trans- 
formations 0{s). U{1) has an inherently geometrical na- 
ture: it is independent of the time rate and only a func- 
tion of the structure of the underlying Hilbert space. The 
non-trivial value of [/(I) is in fact a manifestation of the 
curvature of the Hilbert space experienced by p when it 
is dragged along in the subspace K.. The connection A 
behaves as a proper gauge potential: under a change of 
basis in /C, 0{s) 0'{s) = 0{s)n{s), where [H, rj] = 0, 
A transforms as A^ A' = ^~^AVL + fl"^ Vfl. By choos- 
ing a convenient gauge for which HO' (1)11 = 11, it is 
possible to show that the net effect on the state p ~ pD f 
after a complete cyclic evolution is given by: 



U'{1) = -Pexp 



ds ' 



n 



ds, 



(14) 



which is independent of the gauge chosen Equa- 
tion (|12|) implies that the largest relevant non-coherent 
contribution to the dynamics of Pdf can originate only 
at a further order in the adiabatic parameter rj, i.e. from 
£-1. By substituting equation in and by mak- 
ing again use of the fact that TkpDF = 0, it is possible 
to show that the first order correction to an evolution 
restricted to IC can be given in the following closed form: 



J^i[pdf] = -i[Z,PDF]- 



^{KlkkPDF + PDpAlAk 



(15) 

2AkPDFAl}, 



where = rfe^i = T kD-^G oS and Z = sIrSi = 
iG^g (^D~^^ — D~^) Goff. The super-operator Ci gener- 
ates an incoherent evolution with a typical time scale on 
the order of 77^^. It is worth stressing that, although 
eq. H15|) looks similar to a master equation in the Lind- 
blad form, it generates an evolution that, when restricted 
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to the subspace K., can be non-trace preserving. This is, 
in fact, the case when Hj^A/jII ^ 0. As, by construc- 
tion, the evolution in the whole Hilbert space is trace- 
preserving, this clearly indicates a leakage of population 
from the DFS to the orthogonal subspace, with a time 
rate on the order "q^^ . 

Remarkably, since £i is the largest order incoherent 
term affecting the subspace K. and it is proportional to 
7~^, it implies a decoherence time in K which is inversely 
proportional to the original decoherence time affecting 
the complementary space: a stronger decoherence in /C_l 
implies weaker environmental effects in K,. This appar- 
ently counter-intuitive phenomenon is the essence of the 
adiabatic approximation. It is the fast dynamics associ- 
ated with the incoherent processes acting on K.± which is 
responsible for decoupling the subspace K, from its com- 
plement. This fast evolution averages out the effect of 
Goff and results in a decoherence time scale in IC which is 
quadratic in T . This is the reason why, in spite of the rela- 
tively long time scale T needed for the adiabatic approx- 
imation, the system in K, is guaranteed to be coherent 
for a time which scales with a further order of magnitude 
in T. Alternatively, the same result can be interpreted 
from the perspective of a quantum Zeno effect. The de- 
coupling between K, and K,±^ can be regarded as the effect 
of a continuous measurement process performed by the 
environment. The latter continuously monitors whether 
the system state leaks to and projects it back into 
K, \n. a. time scale of 7~^, which is much faster than the 
leakage process induced by Gog. 

On the light of the previous discussion on £1, a remark 
on the generator of the holonomy is in order. Notice 
that equation H14|l implies that, after a closed motion, 
a non-trivial net evolution is achieved only if i[^,n] = 
Goff + G^jj is non- vanishing. Indeed [^,11] = implies 
dW/ds = 0, which yields to a trivial motion of the sub- 
space /C, and the resulting holonomy has to be trivial. 
This is in agreement with the idea that the holonomy 
is the result of the curvature of the Hilbert space ex- 
perienced by the subspace JC{t) as the latter traverses 
a non-trivial closed motion. Therefore, for an holonomic 
evolution in JC, the latter subspace needs to be coupled to 
via Goff, which, as seen in eq. H15() . exposes po to en- 
vironmental effects. This implies that decoherence on K, 
is in principle unavoidable, as far as a non-trivial holon- 
omy is concerned. On the other hand, as seen before, Ci 
appears as a small perturbation in the adiabatic expan- 
sion, and therefore its effects are adiabatically eliminated 
in the limit of T large enough. 

The results presented in this manuscript reinforce the 
idea that geometric evolutions not only do not depend 
on the dynamical procedures used to generate them, but 
can even be obtained in the presence of strong coupling 
to external reservoirs (non-unitary evolutions). As shown 



here, the interaction with the environment can be used 
to adiabatically manipulate quantum states in a coher- 
ent way, and generate holonomies in arbitrary dimen- 
sional DFS. Apart from their fundamental interest, these 
results, combined with techniques for reservoir engineer- 
ing, may provide alternative approach for robust quan- 
tum computation. As for any other geometrical proce- 
dure, this holonomic evolution has an intrinsically fault 
tolerant nature 13], and due its very nature it is free 
from dynamical contributions, which are known to be 
generally more sensitive to errors. 
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